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1. Introduction 

The super symmetric Yang-Mills (SYM) theories have rich physical content and their 
quantitative analysis is in general as difficult as in the usual, nonsupersymmetric, gauge 
theories. However they often contain, unlike the purely bosonic YM theories, specific 
sectors, which can be analysed exactly and where the supersymmetry leads to a nilpotent 
(topological) symmetry 0. The dimensionally reduced versions of the SYM theory even 
allow various massive deformations conserving this symmetry. 

In refs. |§ and || this symmetry (in the zero dimensional reductions of SYM) was 
applied for the calculation of the (bulk part of the) Witten index for ensembles of N 0- 
branes in 4, 6 and 10 dimensions, justifying the conjectures related to the existence of 
bound states of zero-branes {§ . In [gj the method was applied to study certain correlators 
of BPS states or, in other words, of perturbations of the original reduced SYM theories, 
which preserve part of the supersymmetry. In the case of the zero dimensional reduction 
of M = 1 SYM theory, the large iV limit was studied exactly using the method of @ or by 
the corresponding integrability properties allowing to write an explicit (KP) differential 
equation for the partition function. One of the unexpected results was that, in the large 
iV limit, the physical quantities exhibit an essential singularity at A = 0, where A is the 
coupling of the massive perturbation. The large iV limit of the dimensionally reduced SYM 
theories is also interesting because it may reveal part of the structure of the nonreduced 
theories, due to the Eguchi-Kawai mechanism. 

In this paper we study the one dimensional reduction of the Af = 1 SU (N) pure Yang- 
Mills theory. Unlike the well known and widely used Id reduction to the usual physical 
time 0, J7|, or [RJ, we will retain the "time" along the light cone direction compactified 
on a circle of radius (3. Only this reduction allows the direct use of Witten's localization 
principle PJEl. In order to get rid of the zero modes of the bosonic fields we will deform 
the theory by a massive perturbation corresponding to the 0(2) twisting of the boundary 
conditions on the time circle with respect to a subgroup of the euclidean symmetry 0(4). 
The SYM theory reduced in this way appears to be identical to the compactified hermitian 
matrix oscillator with SU(N) -twisted boundary conditions. The twisting angles are related 
to a global mode of the (time-like) gauge field. This model can be further reduced to that 
of a unitary (already time-independent) twist matrix. We find that the model is integrable 
in the sense that its partition function is a tau-function of Toda hierarchy, i.e. it obeys a 
chain of nonlinear Toda equations. 

The model can be solved exactly and rather explicitly in the large iV limit. The 
solution of the corresponding saddle point equation and its physical consequences in the 
limit of vanishing perturbation represent the main result of this paper. The solution is 
parametrized in terms of elliptic functions. The analysis of the solution as a function of 
the two parameters (3e and A, where (3 is the compactification length, e is the strength of 
the massive perturbation, and A is the twist coupling reveals the following phenomena: 

1. In the double limit e — > 0, A — > the free energy is universal (under certain deformations) 
function of the ratio e/A. 

2. In the limit of vanishing massive perturbation (e — > 0), the observables exhibit an 
essential singularity ~ exp f — co " st ). 



1 Our argument follows essentially the construction proposed in PifLOl], which allows to lift by 
one the dimension of the spacetime without loosing the supersymmetry. Technically speaking, our 
procedure of dimensional reduction replaces the Euclidean spacetime by a point and at the same 
time introduces the "time" dimension r. The latter might be interpreted as a lightlike dimension 
of the original spacetime, but we do not know to what extent this interpretation is justified. 
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3. The analytical continuation of the model at the point e/3 = m (inverted oscillator) 
shows the scaling of the c = 1 compactified noncritical string theory. This value of the 
compactification length e/3 corresponds to the Kosterlitz-Thouless critical point (see ap- 
pendix C). 

The paper is organized as follows. In section 2 we describe the reduction of the 
partition function of the one-dimensionally reduced M = 1 SYM4 theory to that of the 
reduced twisted matrix oscillator, by the use of the supersymmetry and the localization 
theorem. We then reduce the configuration space of the model to the set of the eigenvalues 
of the unitary twist matrix. In section 3 we find that the partition function of our model is a 
r function of the Toda integrable hierarchy and write the differential equations satisfied by 
the partition grand canonical function. In section 4 we give an exact solution of the saddle 
point equation for the large N limit of the model in terms of elliptic parametrization; the 
calculations are presented in Appendix A. In section 5 we study the limit of small massive 
perturbation. We find an universal expression for the free energy in presence of a source 
for the Wilson loops, in the scaling limit e — > and A — > 0. We analyse the properties 
of the solution, especially in the small compactification radius limit and near the curve of 
the Gross- Witten type transitions. Section 6 is devoted to conclusions. In Appendix C 
we give the solution of the analytic continuation of our model to imaginary time e/3 = Z7r, 
which is the Kosterlitz-Thouless point for the corresponding c = 1 noncritical string. 



2. Definition of the model and its reduction to one-dimensional matrix quan- 
tum mechanics 

In this section we will show that the dimensionally reduced M = 1 super Yang- 
Mills theory with gauge group SU(N) can be mapped to one-dimensional matrix quantum 
mechanics. The dimensional reduction consists in replacing the 4-dimensional (Euclidean) 
space-time by a single lightlike " time" . 

Let us first give the generalization of the argument of || to the case of one dimension. 
We start with the the M = 1 SU(N) SYM4 containing 4 bosonic matrix fields (/i = 
0, 1, 2, 3), and 4 real fermionic fields \I/ Q (a = 1, . . . , 4). After performing a Wick rotation 
xq = — 1x4 , the action of the Euclidean theory can be written as 



S = J d 4 xTr {-\Fl v + ~ * T (V 4 + 7 ' V)tt) 



f2.1) 



where V, 



id^ + A^ is the covariant derivative and the gamma-matrices are represented 



as direct products of Pauli matrices: 7$ = Oi x oi (i = 1,2,3). The gauge group SU(N) 
acts to all fields in the adjoint representation, let us assume that all fields depend only on 
the time-like coordinate 

T = £3 — zx 4 , (2-2) 

which parametrizes a circle with radius (3. The resulting model is a matrix quantum 
mechanics containing four bosonic and four fermionic matrix variables. 

We will evaluate the functional integral for this one- dimensional matrix model by 
mapping it to a cohomo logical field theory, which will allow to apply Witten' s localization 
argument [jll[. Let us redefine the fields as 



/ A x \ 








A 2 


x 2 






A 3 + iA 4 
\A 3 - iA^J 




V 


V 4> J 







(2.3) 
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Then the action (|2.1| ) can be written as a BRST exact form. The BRST transformation Q 
acts on the complex of fields $ = {X a , 0, 0, if; ^> a , 77, x} (where H = i[Xi, X2] is considered 
as an auxilliary field ) as 



QX a = i> a , Qip a = [id T + 0, X a ] (a = 1, 2) 
Q4> = 77, Qt] = [id T + 0, 0] 
QX = #, Qtf = K+0, X ] 
Q0 = O. 



(2.4) 



Namely, the action 



S = J drTt(iH[X u X 2 ] + \H 2 + [X a , id T + </>][X a , 0] + \[id r + 0, 0] 2 
+ - X e ab [X a ip b ] + [X a , v ]i} a - [faX a ][id T + (p,X a ] + - X [id T + <j),x] + bPa,4>}ipa 



can be written as 



rP 

S =Q rfrTrV($) 

2 (2-5) 

v($) = i v [td T + 0, 0] + x (h -i[x 1 ,x 2 }) + J2 M*a, 4>}- 



a=l 



The square Q 2 of this transformation represents the "time" gauge transformation generated 
by 0, V r = [id T + 0, ] . Hence Q is nilpotent on the gauge-invariant quantities. The ghost 
number of the fields is —2 for 0, —1 for 77 and x ? for X a and if, +1 for ift a , and +2 for 0. 
The functional integral with respect to the BRST complex of fields $ 

(where the integration measure is normalized by the volume of the gauge group Q) can be 
therefore evaluated using the Witten's localization argument ||11|| . Namely, the integral is 
saturated by the BRST critical points Q<& = 0. More strictly, we have to integrate over 
a continuous critical manifold, because of the zero modes of 0, X , 0. The zero modes are 



elliminated by adding, following a Q-exact term to the action ( |2.5|) by changing the 
action to S + SS, with 

fP 

5S = tQ / rfrTrx0. (2.7) 
Jo 

We can also discard from the very beginning, the term ^rj[id T + 0,0]. As before, H can 
be integrated out by setting H = t<p + [Xi,X2] in (|2.5|) . The advantage of introducing 



the perturbation ( p.7|) is that for t ^ 0, the fields 0, X , and 77 can be integrated out. 
However, the perturbed integral does not coinside in general with the original one because 
of the new fixed points "flowing in from the infinity " when one perturbes to t 7^ 
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[ 1 1]. The correct statement is that there exists a class of BRST-invariant operators whose 
(nonnormalized) expectation values coincide in the original and the perturbed theory. The 
first such operator is 

rP rP 
u = / drTr(-[id T + 0,Xx]X 2 + fafo) = Q / dre ah lh:^ a X h ). (2.8) 
Jo Jo 

The second is any 5'C/(iV)-invariant function /(fl) of the holonomy factor around the circle 

O = Texp (i [ P dT(p(r) ) . (2.9) 




The functional integral for the expectation value 

{eWfm = I WQ) e ~" iS+5S)+LUf{n) (2 ' 10) 

does not depend on the coupling t. Indeed, taking the derivative in t and integrating by 
parts, we find zero, since the integrand vanishes at infin ity d ue to the factor e u . Therefore 
we can take the limit t — ► oo, after which the integral ( P-10|) gets localized near the zeros 
of H, cj), x, r). (In particular, the partition function does not depend on the gauge coupling 
g.) Now we can calculate the partition function ( |2.6| ) understood as the average of the 
identity operator. Since u has ghost number +2, due to the ghost number conservation 

Z N (/3,g) = (l) = (e u ). (2.11) 

The above argument has been applied recently by F. Sugino [|12] in order to calculate the 
partition function of the four-dimensional M = 1 SYM reduced to a two-dimensional torus, 
with periodic boundary conditions for all fields. Our case is slightly more subtle, because 
of the zero modes of the fields X a . These zero modes will be elliminated, as in [Q and 
later in |J, namely by deforming the BRST operator in the definition of the action ( |2.5| ). 
Let us first notice that after the redefinition of the fields the theory is still invariant under 
the 0(2) rotations in the directions orthogonal to the light cone: 

Xx + iX 2 -> e ie (Xi + iX 2 ) , tpi + itp 2 -> e ie (^i + # 2 ) ■ 

This allows to construct another BRST operator, which squares to a linear combination 
of a gauge transformation and an 0(2) rotation. The twisted BRST charge Q e acts as 



Q e X a = ip a , Qeipa = [id T + (/>, X a ] + iee^Xp, 
Q e cp = 77, Q e r) = [id T + 0, 0] 
Q eX = H, Q e H = [id T + cf),x] 

Qe4> = 0. 



(2.12) 



The modification of the supercharge is equivalent to changing the action ( |2.5|) and the 
operator (|2.8|) as 



S -> S + 2ie J dTTr(<t>[X u X 2 ]) 

/ drTv(X 2 1+ Xi). 
Jo 



e 

u — » u 

2 



(2.13) 
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In the limit t — > oo the integral gets localized near the zeros of H, <f>, x, ?7, leaving the place 
to the action 

S = J drTr (-i[id T + fa X X ]X 2 - \e{Xl + Xf ) + ipiip 2 ) , (2.14) 

and the ^'s can then be integrated out. Finally, the integration over X 2 gives the partition 
function of the matrix oscillator (with the coordinate X\ = X) in presence of the one- 
dimensional gauge field far) 

Zn(P, 9,e) = J m y^g (T) exp (~Tr j" dr Q [zd r + fa X? + eX^ j , (2.15) 

with periodic boundary conditions X((3) = X(0), fafi) = 0(0). It is clear that the integral 
depends on (3 and e only through the product e(3. We will absorb e in /3, 

ef3^(3, (2.16) 

remembering that the perturbation is lifted in the limit f3 — > 0. 

The functional integral over the field (p can be written, after fixing a gauge d T (p = 0, as 
an integral over the unitary matrixH representing holonomy factor defined by ( |2.9| ) , namely 
O = e l/3< ^ normalized by the volume of U(N). The holonomy factor enters the functional 
integral over X as the twisted boundary condition: 

X{J3) = tt+X(0)tt. (2.17) 

The integral over X can be performed exactly, and the integral over the unitary matrix Q 
reduces to an integral over its eigenvalues e %Ql , e l0N (which are defined up to a permu- 
tation, hence a combinatorial factor 1/N\). The partition function is therefore given by 
the iV-fold integral! 

where 9i = (3 fa. 

The expectation value we are calculating is a deformation of the Witten index 

Z N (J3) = Tr(-) F e- /3H e l/3eJ , (2.19) 



2 It is assumed that the integration contour for the eigenvalues of <f> is chosen along the real 



axis. In this case <p should be taken anti-hermitian, see the the discussion in [11|. 

3 This happens to be exactly the partition function of the one-dimensional gas studied by 
Michel Gaudin in 1966 JTl; it was extensively used in |L4j to study the compactified 1+1 dimen- 
sional string theory via matrix quantum mechanics; in relation to the actual SYM theory this 
formula was communicated to us by N. Nekrasov. 
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where the (— ) F -factor is included in order to impose the periodic boundary conditions on 
the fermionic fields, and the trace is twisted by an 0{2) rotation e l @ eJ in the (12) plane. 
The twisting of the BRST charge Q — > Q e does not change locally the functional integral, 
but it does change the boundary conditions for the fields. 

Now we are at the most subtle point of the reduction procedure, which deserves to 
be discussed in more detail. Considering the fields <f> and 4> as two independent fields 
imply that the integration over them is understood as contour integration. The twisting 
separates the poles and zeroes of the integrand and allows to evaluate the integral by the 
residue theorem, in complete similarity with the calculation of for the zero dimensional 
model. Since the integrand does not depend on the variable 

- 9 1 +... + 9 N 

°= N ' 

the contour integral with respect to this variable would give zero. In fact, the integration 
with respect to this variable should be excluded because this is one of the normalizable 
zero modes of the original fields and the measure should contain a product of delta 
functions of the bosonic and fermionic zero modes. In particular, the normalized zero mode 
of is 

,/3 



For a more detailed discussion see [jI5| . Therefore the measure in ( |2.18| ) contains a delta 
function 5(<p^) ~ 5(9i + ... + On), which suppresses the contour integration with respect 
to 9. The integral over 0's is normalized by the volume of the residual global gauge group. 
The introduction of the delta function should respect this normalization. Thus we have to 
insert 

2n5(9) = 2nN5(9 1 + ... + 9 N ). 

Now we can integrate, after representing the integrand as a determinant using the Cauchy 
identity, by using the residues theorem. The integral is equal to the sum of the identical 
contributions of the (N — 1)1 cyclic permutations in the expansion of the determinant 

ZnV) = R^^iH / ^...^2nN5(9 1 + ...9 N )f] — T7 l - - 

K > y> N\ J 2n 2n 11 NJ 1 = 1 S in[§(^ - 0<_i + ift] 



I 



dOi d,0 N fi/Aj Oi+ON 



sinh \[9x -9 N - i(N - 1)(3] sinh \[9 X -9 N + i/3] 
1 



2A^sinhiVf 



(2.20) 

In the limit (3 —>■ oo our partition vanishes, which is not unexpected, since the Witten 
index of the D = 4 theory is zero. In the limit j3 —>■ we recover the result for the 
completely reduced theory ~ 1/iV 2 , in agreement with [Q. A more careful analysis allows 
to reproduce also the numerical coefficient, in accordance with the conjecture made in (Til 



In the limit j3 — > 0, the ^-integral is saturated by the integration in the vicinity of the 
iV central elements of SU(N) which are parametrized by the element of the su(N) Lie 



algebra [I7|. After performing carefully the limit, one finds (see e.g. 
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where zff(g) is the partition function of the completely reduced theory, and one repro- 
duces the result of 0: 

zf{g)=F N9 -^ 2 -vl- 2 . (2.21) 

(The numerical factor Tn depends on the way the integration measure is normalized. In 
the normalization used in || this factor is equal to one, but this is not the most natural 
choice from the point of view of applications to the D-brane physics.) In the particular 
case of the SU(2) theory the results was obtaind by the direct calculation of the integrals 

Let us note that our partition function only formally coincides with that of the twisted 
matrix oscillator and, at least for finite N, there is an ambiguity related to the prescription 
for the contour integration. Witten's localization procedure leads to an integral over the 
Lie algebra and logically the integration with respect to 0's should be taken along the whole 
real axis. With this definition, only (N — 1)! terms in the expansion of the determinant 
will contribute to it. On the other hand, had we integrate in interval [0,27r], this would 
correspond to contour integration with respect to the eigenvalues of the Lie group element 
{tk = e' l9k }k = i, where the contours circle the origin. In this case we would get contributions 
from all AH terms in the expansion of the Cauchy determinant. The result would be given, 
instead of ( |2.2U|) , by (see, for example, ref. JH[ ) 



e -N 2 p/2 



Zn(P) ~ (l-e-^)(l-e-^)...(l-e-^)- (2 ' 22) 



Unlike (|2.20|) , the (5 — > limit of this formula does not match the result of 



Which of the two formulas ( |2.20| ) or (|2.22[) is correct? Clearly the difference between 
t hem is due to a different treatment of the boundary conditions for the field <p in the formula 
( |2.15|) . result considered then at hermitean oscillators). A happy resolution of this paradox 
would be that from the point of view of the application of Witten's localisation principle 
both formulas seem to be possible but the result depends on the bo und ary conditions and 
the contours of integration for the field <p(t) in the original action ( |2.5|) . However we feel 
that the question is rather subtle and more study is needed to clarify it. For example we 
cannot be sure that the supersymmetry of the original model isn't violated in one of two 
cases. On the other hand, the local BRST symmetry used for the calculations is certainly 
intact. 

Now let us consider a slightly more ambitious problem, namely to calculate the gen- 
erating functional of a set of BRST invariant operators made out of the gauge field 4>. 
As mentioned before, such operators can be constructed as traces of the holonomy O in 
different representations, or, equivalently, as polynomials of the moments Tr(0) fc . We will 
add to the action the simplest possible source term 

ATr(0+ + 0). (2.23) 

Repeating the arguments, which led to (|2.18| ), we find for the generating functional the 
following integral representation 

N] J fc=i 2n tfj Sm k \ Q i ~ 9 J + 

where 6i = (3<pi. If one follows the recipe of ||, the integration should be considered as 
a contour integration along the real axis, where the A" integration variables are subjected 
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to the constraint 9± + ... + 9n = 0. Then the result should be analytic as a function 
of P, which can therefore be given complex values. It is plausible that in the large N 
limit, which we are interested in, if the perturbation is sufficiently strong, the choice of the 
contours should be not important. The equivalence between the reduces SYM theory and 
the twisted matrix oscillator should takes place only in this limit. 

It would be very interesting to understand what is the meaning, in terms of the original 
supersymmetric theory ( |2.1| ), of the deformation that leads to the partition function ( |2.23| ). 

The reduction from 4 to 1 dimension of the original theory turns three of the compo- 
nents of the gauge field into Higgs fields (Xi,X 2 and <fi). This makes the direct calculation 
of the partition function (which is related to the bulk part of the Witten index) more 
delicate, because of the absence of mass gap. By introducing the deformations (|2.7| ), ( |2~8| ) 
and (|2.13| ) we add an additional Higgs potential, thus breaking part of the supersymmetry. 
The effect of the source term, which we added to obtain the partition function ( |2.23| ) , de- 
pends substantially on the way we have perturbed the th eory Indeed, it has positive ghost 
charge, and its effect would be zero, if the perturbation ( |2.13|) of ghost charge —2 were not 
there to compensate it. This is also true in the completely reduced theory, discussed in 
ref. H. 



3. The partition function as a tau-function of the Toda hierarchy 

Here we will show that our partition function with a source term ATr(0 + + O) is a 
tau -func tion of discrete Toda chain. Let us rewrite the partition function of the model 
eq.( p.24|) in the following form: 

where Zk = e lBh , U(z) = J^n^o^nZ 71 > an< ^ ^( z ) * s the Van-der-Monde determinant of z's. 
In our case t\ = t-i = iVA and t n = for n/±l, but most of the following conclusions 
are true for a general U(z). 

Let us now introduce the grand canonical partition function with the "charge" /: 

oc 

n[t,fi} = e» N e- lN ?Z N {M (3.2) 

N=\ 

Due to the Cauchy identity the last equation can be rewritten in terms of a functional 
Fredholm determinant: 

*M=xy" v/n^ as ^ (3 .3) 

= Det(l + e»-P l K), 
where the operator K is defined as 

r dz e wt z )+u{z')] 
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It is convenient to modify slightly the definition of the tau-function: 

Ti[T,fj] = fi[t, fi]exp(- J^nt n t_ n ), (3.4) 

n>0 

where we introduced new couplings T n by: 

U(z) = z n T n {e- np/2 - e n/3/2 ) (3.5) 

so that the old couplings are expressed through the new ones as: 

t n =T n (e- n ^ 2 -e^/ 2 ). (3.6) 

We also note that 

n[T, fjt] = t [T, n-pi] = r[T, n-/3l\. (3.7) 



Using, for example, the general construction of the paper []19j for our particular tau- 
function we conclude that it is a particular case of the tau-function of Toda hierarchy. It 
satisfies the Toda chain equations. Namely let us introduce a new function 



$, 7l-i r[T, fj, - p{l - 1)] 



n t[T,/j, — (31] [ - 



and the notations 

ix± = V2T ±1 = ±(e^ 2 - e-V 2 )-H ±1 . 
The first equation of the Toda hierarchy can be written as 



() 3 + ! ( e *«-*«+i _ e *i-i-*i) =0 (3.9) 



dx+ dx- 2 

Due to the symmetry 9 — > —0 of the measure our tau-function depends only on the 
variable x = y/x+x-. The corresponding reduced equation is: 

$" + + ! ( e *i-*t+i - e *i-i-*i) = o, (3.10) 
x 2 

where the derivatives are taken with respect to 

For the function ^ = — = log r '~ 1 I' +1 the Toda equation reads: 

¥i + ~i>\ + ~ (2e^ - e^" 1 - = (3.11) 

x 2 

The tau-function, as well as <&z(0) and -0^(0), can be determined for x = using the 
methods of [irj| (for x = the tau-function is the grand canonical partition function of the 
matrix oscillator in the singlet representation of the U(N) group, which is the same as the 
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partition function of iV fermionic oscillators) and it can serve as a boundary condition for 
the Toda chain equation. For example, one finds from ( 2.22 ) 



e M0) = 1 + 6 (3.12) 

Let us also note that if)i(x) is analytic in x 2 at the origin, which gives the second initial 
condition d x ijji\ x =o = 0. 

Using these equations and the boundary conditions we expand the partition function 
( |2.24|) in powers of A 2 . In the first order: 

1 d 2 , _ 1-e-^ 



2N*d\ 2l ° gZ = T^- (3 - 13) 

This is the simplest correlation function (trfi+trfi) of the holonomy wilson loop in 
our original model. 

The large N limit of the initial partition function (|2.24 ) or (O) can be studied in 



terms of a special scaling limit of these Toda equations (since \i ~ N in the legendre trans- 
form from canonical to microcanonical partition function), similar to the KP-hierarchy 
approach of a simpler zero-dimensional model of paper |j . We leave this study to a future 
publication. 



4. Saddle point equations in the large N limit 

In this section we will investigate the large iV limit, which is the most interesting from 
the point of view of applications. Since the potential A cos 9 is symmetric, we assume that 
the saddle-point spectral density 



1 N 



i=l 



is supported by the symmetric interval [—a, a] with < a < it. The function p{9) is 
determined by the saddle point equation 

2Asin0 = pe'p(e') (2 cot e -^L _ C ot d ~ d ' + ^ _ cot 

where we temporarily rescaled (3t — > (3. This equation is equivalent to a functional equation 
for the resolvent 

W{9) = \ f de'p(O') cot (4.2) 

namely 

\sm9 = W{9 + i0) + W(9-i0)-W(9 + i(3)-W(9-i(3) (4.3) 
where —a < 9 < a, supplied with the normalization condition for the density 

±W{z) = 1 (4.4) 
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where the contour of integration C circles interval [—a, a]. 
It is easier to solve this equation for the function 



—2 cos z — 



A i 
4sinh 2 § f a dOp{d) 



J —a 



(4.5) 



A J_ a cos(z — 9) — cosh | 
which satisfies the simpler equation 

C(* + if) = C(*-i§) (»e[-a,a]). (4.6) 

The solution can be formulated in terms of standard elliptic functions (see Appendix 
A for the derivation). We give it in the form which is convenient for the limit of small (3 
(or, equivalently, finite (3 and e — > 0). The function ((z) will be given in a parametric form 

C = ((v), z = z(v) 

where the parameter v belongs to the rectangle — § < Ret> < f , — § t < Imv < ^r. The 
elliptic modulus q and the nome k 2 

q = e -*K/K' = 

7 _ yj f l-q^ Y (4-7) 

k ~ 11 ^l + (? 2n-lj 

n=l x J ' 

are given below as functions of (3 and A. 
The solution (in parametric form) is: 

_ C4/ 2 (^)-C3/ 2 (^) , 4 ^ 

/ 2 (v)-f 2 (v oc ) 

where f(v) is a standard elliptic function 

/(«) = — dn(— V , fc'j =1 + 4^^-^008(2^), (4.9) 

n=l 



and 



(4.10) 



= i-i> + 1 In — r 

./5 sin(w + v ) , . o 2n sin 2m; sin2nw 

= i-v + iln ) 2£Z + 4 Z \ * 22 . 

7T sin(i> — v ) ^-^ 1 — q Zn n 

The modulus r is proportional to the ratio of v x and (3 

r = 4z^, g = e~^ t, <», (4.11) 
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and is determined by 



2 i sinh(/?/2) = lE(k) - «£±<l£iif (*). (4.12) 



The parameters (j, £5 of the solution are expressed as functions of A and v x as follows: 



/(«-) »i(0) 



e^oo ff"( v ) f'( v ) , fl 



(4.13) 



Cl " «2-Jfe2 ' Cs " l-a 2 ' (4 ' 14) 

a= f ~M „ 7= ; (C4 " C3)Q (4-15) 

l + 2E"if v/(l-a 2 )(a 2 -P) 1 J 

Finally, it is useful to know the value (2 of the function £(z) at the branch point z = i^ +a 



C4-C2 P , P ^ • o , lfil 



7 27T 

' 00 n=l 



5. Scaling limit 

5.1. The resolvent in the scaling limit 

Let us recall that the parameter f3 is the product of the physical time and the twisting 
parameter e. Therefore the twisting is removed in the limit /3 — ► 0. If /? — > with A fixed, 
we reproduce the zero-dimensional case considered in ||. In this section we will consider 
a nontrivial limit where both A and f3 go to zero so that the ratio X//3 remains finite. In 
this limit, all observables depend only on the ratio \//3, and this is why we call it "the 
scaling limit". Note that in the thermodynamical limit N — > 00, the two limits A — > and 
e — > will not commute. 

In the scaling limit we have Imr >> 1, K^ln^j, K' = §(1 + ^r) and, neglecting 
the exponentially small terms, we get 

, = C4 + C3 C43 1 



where (Q k = (i - Ck) 



2 16q sin(u + O sin ( v - O , g ^ 

7T sm(f — f oo ) 



(43 = 16qe^ v °° sin 2?;^ 

C4 + C3 t v f0 L /? . ( 5 - 2 ) 
= — 00 (2 cos 2i> H sm2t> I. 

O V 00 1 00 / 

2 7T 
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When the regularization is removed, i.e. in the limit (3 — > 0, a sensible limit is obtained 
when A tends to zero linearly with (3. The scaling coupling constant (3/X is obtained from 
((TLJ) after substituting E = 1, K = ^-v x : 

27r^ = 2sin2w - 4v cos2i> . (5.3) 

5.2. The free energy in the scaling limit 
The derivative of the free energy 

F(A,/3) = hm^lnZ N (A,/3) 

is proportional to the first moment of the spectral density 

F' x (\,(3) = ! d9p(9) cos 6, (5.4) 

J —a 

which can be evaluated by looking at expansion of ((z = n + iy), at y — > oo, 

oo 

C(z) = e» + j^C (fc) e-^ k +^. 

k=l 

We have 

F' X (\P) = — ^ [C (1) -l]- (5.5) 
8 sinn ^ 

The coefficient is evaluated in the scaling limit (3 — > in Appendix A. This allows us 
to write an explicit expression for the free energy in the limit f3 — > 0: 

F' x (X/[3) = ^-(4v ac -sm4vJ, jrf = sm2 Voc - 2v x cos 2v x . (5.6) 
47rp A 

This expr ession for the free energy is universal in a certain sense: if one deforms the 
potential (|2.23|) to a more general one: 



ATr J2 ^° n ' ( 5 - 7 ) 

n= — oo 

then the scaling limit of the free energy will have the same form ( p.6| ), where A will be 

substituted by some function of the couplings X(gi,g2, ■ • •)• The universal form of the free 
energy can only change if we tune the couplings g n to some multicritical point. 

The corrections to the eq. (|5.6| ) are of two kinds: power-like corrections and exponen- 
tially small terms of the type 

( 4X 1 / 3 

In the limit A — > oo we have q = exp — ( J . These terms are of course invisible com- 
pared with the power-like corrections but they imply the existence of essential singularity 
in the (3 -> 0. 

If we return to the original notations in terms of /3,e and A we conclude from (5.6) that 
F(X, f3, e) = eXf(X/ {(3e)). Hence the principal ~ A^ 2 correction to the free energy tends to 
zero in the limit e — > (when we recover the original unperturbed reduced SYM theory). 
On the other hand, as we will see below from (5.6), there is no regular expansion in powers 
of e in the weak coupling phase which signifies that there is an essential singularity at the 
origin of this coupling and, correspondingly, in the moduli space of our theory. 
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5.3. The Gross- Witten phase transition and the strong coupling phase 

The matrix integral we are considering has qualitatively the same phase structure as 
the U(oo) gauge theory on a two-dimensional sphere. The weak coupling phase considered 
above, describes the range of couplings A > A c where 

1 -e~P 

X c = (5.8) 

is determined by from the condition a = n, i.e. that the two endpoints of the cut meet on 
the unit circle. (Eq. ( |5.8|) follows from (|5.3|) with = 7r/2; the length of the cut as a 
function of A is given in Appendix A.) The singularity near this critical point is as usual 
of third order. 

The strong coupling solution is obtained by expanding the spectral density and the 
kernel in a Fourier series 



p{9) = — } c k cos(k0). 

n=0 

6 + iO 6-iO 9 + ifi 6-if3 



cot — ; h cot — cot — — - - cot — — - =4V[1- e~ kf3 } sin(fc0) 

£ £ £ 



k=i 



de'p(e') cot = J2 Ck sin ( ke ) 



k=i 



It is therefore clear that only the Co term of the expansion of the spectral density has 
to be retained. One finds 

p(0) = -L(l + A _. cosfl) for 0< A< A c = 1 - e _/3 . (5.9) 
Z7T \ 1 — e " J 

For the free energy we find then: 

F' x = f d9p{9) cos# = | _ A _ for < A < A c (5.10) 

J — IT ^ 6 

In the scaling limit e — > we obtain: F' x = |^ for A < (3. At the critical point A c = (3 
we have F' x \\ = p = | and F x \\ = p = A simple calculation using eq. ( |5.6| ) gives in the 

weak coupling phase A > A c the same values of first two derivatives of the free energy at 
the critical point. This means that we have, as usually, the 3-rd order Gross- Witten phase 
transition. Note that in the limit (3 — ► oo our model reduces indeed to the one-plaquette 
model originally studied by Gross and Witten [pQ| . 



5.4- Reduction to the zero- dimensional theory: (3 << A 

In this limit the theory appears to be the zero-dimensional reduction of CN = 1 SYM 
studied in . The integral ( |2.24j) reduces (after the rescaling 9 — > (39) to a simpler integral: 



k—l i^j 
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This model was studied in || and later in 0. 

We find from (|5.6| ) the following expansions in half-length of the cut a = 2v c 



2 4 

er a 
~ 10 + 4200 



+ 0(a 6 ) 



a 



A 3tt x " 10 
This gives the following asymptotics for the free energy: 

2/3 



Ftf, A) = A 



3 /37T/3 
10 



A 



+ 0(a 7 ) 



27 /3tt/3 



1400 



A 



4/3 



+ o 



(5.12) 
(5.13) 

(5.14) 



The first term of this expansion matches the asymptotics of big A obtained in || for the 
integral Q5.ll ) in the large iV limit. The next terms are not supposed to match with 
|| since we already used the scaling limit expression of the free energy with the finite 
compactification radius (3. 



6. Conclusions 

Let us outline the main results of the paper: 
1. We consider a topological sector of the one-dimensionally reduced Af = 1 SYM4 theory 
on the light-cone time circle with a a special massive perturbation. Using Wittens non- 



abelian localisation principle [11]], we represented the partition function with periodic b.c. 
in terms of a solvable matrix quantum mechanics (twisted matrix oscillator). 

2. We find the integr ability properties of this model relating it to the To da hierarchy. The 
generating functional as a function of its parameters satisfies the Toda chain equation. 

3. In the large iV limit we find the exact solution of the model: the generating functional is 
parametrized in terms of elliptic functions. We find the Gross- Witten type phase transition 
and identify its location. The strong coupling solution is also found. 

4. An interesting model corresponds to the analytical continuation /3 — > i/3 (inverted 
matrix oscillator) being known to have the properties of the c = 1 non-critical strings. 

5. In the scaling limit of vanishing perturbation we find a simple universal (with respect 
to certain deformations of parameters of the generating functional) expression for the free 
energy and Wilson loop correlators along the light-cone circle. Its strong compactification 
limit restores similar results for the completely reduced M = 1 SY M4 considered in || . 

Some remaining problems: 

1. We need further understanding of the space time symmetries of the model and of the 
correlators corresponding to our generating functional. 

2. The representations similar to the eq.( |2.24| ) for the SYM4 can be found also for the 
SYMq and SYM\q reduced to the light-cone time circle: we just have to take the cor- 
responding eigenvalue integrals for the partition functions in the paper [||and substitute 
there the rational functions by trigonometric ones. Unfortunately, we cannot apply the 
powerful methods used here to those models: we don't know any relation of them to the 
integrable hierarchies and we cannot solve exactly the large iV saddle point equation. On 
the other hand, it seems to be possible to investigate this saddle point equation in the 
scaling limit similar to tha used in the present paper. 

Another interesting question is whether two prescriptions for the contour integration 
with respect to the eigenvalues of cf> coincide for the infinite N in some part of the phase 
space of parameters e(3 and A. It is clear that they give different answers in the strong 
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coupling phase A < A c (/3) (since they are already different for A = 0). As for the weak 
coupling phase A > A c (/3) of our large N solution (which is not even analytical at A = 0), it 
is possible that the saddle point approximation does not distinguish between two different 
prescriptions of integration over cj)(t) (contour integration over the Cartan subalgebra, on 
the one hand, and integration over #'s in the finite interval [0, 2tt], on the other hand). 
This hypothesis is to be verified. A w eak er version of it could be the coincidence of two 
prescriptions in the scaling limit (eq. ( ^6|) ). 
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Appendix A. Solution of the saddle point equations 

A.l. The function z = z(Q as an elliptic integral 

It follows from the integral representation of ((z) that it is real when z G H, iH, zIR±7r, 

satisfies 

C(^) = C(^ + 27r)=C(-^)=C(^) (A.l) 

and by (fO|) is also real along the interval — a, |/3 + a]. Therefore this function defines 
a map of the half strip < Rez < 7r,Im^ > with a cut |/3 + a], to the upper half 
plane Im£ > (Fig. 1). The inverse map z = z(Q is given by the Schwarz-Christoflel 
formula (see, e.g. [plj): 

■'*(*-<»> (A.2) 



C4 



Y(t) 



where 



Y{t) = v / (^-Ci)(t-C3)(t-C 4 )(t-C5) • (A.3) 

By construction, the map (|A.2| ) acts on the special points Ci < C2 < C3 < C4 < Cs an d 00 
as is shown in the two first coloumns of Table 1. 

The values of ( at the special points of the map are determined as functions of (3 and 

A by the assymptotics of £(z) at infinity. The expansion of the function ( |4.5[ ) at z — > 00 
contains only odd powers of e lz . If we approach infinity as z = tt + iy, y — >• 00, the 
asymptotics of £(z) is 

00 

C+(y) = + iy) = e y + J2 C (2n+1) e"^ +1 ^ (A.4) 

n=0 

where 



^(1) 



8 sinh.^ — 

1 + — — / d0p(6) cos6 \, (A.5) 



etc. 
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A. 2. Elliptic parametrization of the solution 

The map ( |A.2| ) and the condition ( |4.4| ) can be expressed explicitly in terms of standard 
elliptic integrals (see [B2[, 256.02 ) with parameters 



C54C31 ,/ n. /C43C51 



fc = i^, k' = Vl-k* = Jff± (A.6) 

S53S41 V S53S41 



7 = VC53C41 , « 2 = 7^, z/ = arcsinW^- (A.7) 

C.53 V S51 



where the notation dj = d — Q is used. Namely 

2C43 ( tt, . 2 7 \ , C32 



7 V C43 



(A.8) 
. 1 C-C 4 

(f = arcsm — \ / — , 

a V C - Cs 

and 

f S inh(/3/2) = 2j5(») - CA + ClC, Jf(t). (A.9) 

It is convenient to introduce as a parameter the elliptic amplitude w related to the angle 

<£> as snu = sin 99 



■ ^ f C dt .1 /C-C4 in v 

«=*^— ] U YW) =FM ' SnU = SmLP= a\j— 3 - (A ' 10) 

Then the function 

C 4 - C3Q 2 sn 2 M 

C u ) = — i — ^—2 — ( A - n ) 

1 — arsn^it 

maps the upper (^-half-plane is mapped to the rectangle < Keu < K, < Im« < K' with 
K and -PC' being the complete elliptic integrals associated with the moduli k and k' . The 
special points C = Ci> • • • > Cs an d 00 correspond to the points u\, . . . , 115 and «oo along the 
boundary of the rectangle as is shown in Table 1. Note that 

- = sn(uoo). 
a 



c 


z 


U 


V 


— 00 


+ioo 


Uoo + i0 


Voo+0 


G 


if* 


K + iK' 


7T 

2 


C2 




u 2 


V2 


Cs 




iK' 


f(l+r) 


C4 










c 5 


7T 


K 





+00 


7T + ZOO 


Uoo + i0 
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Table 1. The values of z, (, u and v at the special points of the map. 



The function z(u) defined by the integral (|A.2| ) reads, in the parametrization ( |A.10|) , 



z{u) 



X 



43 



41 



o 1 — a 2 sn 2 u (43 



(A.12) 



We will express the integral in ( |A.12| ) in terms of Jacobian elliptic functions. Since the 
point is between u$ and ui), it has the form 



K + i£, 0<£<K'. 



We find, using eq. 433.01 of 



z(u) = 1 In -^-^ — + — 



X 



4.2 



H x {u-i£) K[^C^i 
where H\{u) is a standard Jacobian elliptic function and 



A (u, k) = — [(E - K)F{v, k') + KE(u, k')} 
n 



K - 7rA (^, k) 



(A.13) 



(A.14) 



■ H'M) 

'#i(*0 



z/ = arcsm 



'c 



41 



Va 2 - k 2 



arcsin 



is known as the Heuman's Lambda function. The condition z(iK') = z(K+iK') is satisfied 
only if the coefficient in front of the linear term in u is zero, hence the condition 



K 



42 



V C53C41 

From H\(u + iK') = e~ ~ H\(u — iK') we find 



K = 7rA (^, k) 



(A.15) 



7T 



z(iK') = i-£ 



which allows to determine £: 
The final expression for z{u) is therefore 



z{u) = i In 



H 1 (u-i§(3) 



(A.16) 



(A.17) 



A. 3. The dual modulus 

We are going to write our solution in a form, which will allow to perform painlessly 
the scaling limit (3 — > 0. In this limit 043 ~ — i(3^ — > and, according to (|A.6| ), k' ~ 
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Ae K — > 0. Therefore it is more convenient to expand the solution in the dual modular 
parameter 

q = e -* K / K ' = (A.18) 



and use the variable v 



n tt/2 

V =- T - % —U 



as a parameter. The elliptic nome is expressed as a function of q as 



n=l 



1-q 



2n-l 



+ q 



2n-l 



(A.19) 



The parameters corresponding to the special points of the map are given by the last 
coulumn of Table 1. The parameter of infinity is equal, by (|A.13| ) and ( |A.16| ), to 



K lf > = -%p (o<„_<£ 



AK 



(A.20) 



The domains of the four variables z, C, u and v are depicted in Fig. 1. 

We will write the solution as a function of the parameters v x = ^(A) and 3. It will 
be written as a series in the expansion parameter q 



_ 4tt , 

q = e p 



r = Ai- 



3 



(A.21) 



The expansion of the function ((u) is obtained by plugging in eq. QA.11J ) the representation 
of sn in terms of the dual modulus 



1 



snu 



(hi ( v, k' 

TV 



7T 



2K 



7 /(?) 



(A.22) 



"t: — cos(2nt>) 
— ^ 1 + q 2n V ; 



(A.23) 



The function z(v) reads, in terms of the standard elliptic functions associated with the 
dual modulus, 



4 fiUv + v ) 



2n s in2nf sin2m> 

J- 00 



3 , sm(v + v ) 

i-v + iln — ^ 

7T sm(v-i; J t^i l ~ q 



2n 



n 



(A.24) 



Finally, ( |A.15|) expands as 



Av C42 2 ) 2 



,2n 



sin 2m; „ . 



n=l 



(A.25) 
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In order to fix completely the solution, let us consider the vicinity of the point v x and 
compare the explicit dependence ( = ((z) with the asymptotics ( |A.4| )at z — > tv + ioo. The 
half-line 



z = 7T + iy (y > 0) 



is parametrized by the interval < v < v x . In the left vicinity of the point v o 



v = v oc -e (e > 0) 



the functions z = n + iy and £ have the form 



A 



P 



e y = - + B + Ce, (=- + Q + Re 



with 



p = -c 



43 



/0>c 



2/'(n, 



g = -p 2 



q />J W i QO 

'C43 f(vj 2 /( V J 2 /'( V J 



a« 0i( 2 « 



^(0) 



P = —A — + n ° 



0l(2 V J 



The leading asymptotics ( |A.4j ) of £(2) is achieved if A = P and B = Q, which yields 

(EH)- 

^From ( |A.25|) we get 



rz — : — P ( sin 2v 
C 42 = VC^Uitt 1 + 27rg s 

27T V V 



(A.26) 



(Note that the relation 2^ = Ci + C3 + C4 + Cs is satisfied.) 

20 





V 



Fig.l: The domains of the variables z, £, u, v. 

A. 4- The limit of large r (small (3) 

In this limit, which corresponds to the scaling limit discussed in Section 4, the v- 
rectangle can be replaced by an infinite half-strip and elliptic functions degenerate to 
trigonometric functions. After substituting E = 1, K = jfV^ in the normalization condi- 
tion ( |A.9| ), we get in this limit 

-4- sinh — e~^ Vco =2sin2v — Av \cos2v — — sin2t> 1. (A. 27) 

\ o ckj oo l oo o oo J V / 

A z Ztv 



The parameters of the solution are obtained from 

Q = C 4 + Pcot2u ooJ R = P I - + cot 2 2v c 



C43 1 ^ /- 1 7-1 L ci-. r> n f 2 , . /2 



16q sin 2?;^ 
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A = e* v °° sin2?; , B = -A (cot 2u + - V C = A (~ + - cot 2u + 



W ' V 3 ' 7T "°° ' 27T 2 

From i = Pwe get 

C 43 = 16qe« v °° sin 2 2v oo (A.28) 

and B = Q implies 

= -el"- (2 00s 2« +-sin2z; ) = - fe-"- sin2« V . (A.29) 
2 7T V "V^ 

It is useful to note that 

C53 + C54 ,2 C42 + C32 /3 C43 i« 2 
— = cot v , — = — cot v , — = Ibq cos v . 

C41 + CC1 X C41 + C31 27T " C41 + C31 

Finally, the coefficient is obtained as 
(W=A(R-C) 

=1 I if [ 4 ^" S ' n4 ^ 1 " sin22 ^ -2^ sin 4^ | (A.30) 

A. 5. The length of the cut 

The branch point of the Riemann surface of ((z) is at Z2 = zfa), where C( v 2) = (2- 
Taking the limit of ( |AT26|) , 

C42 = -sin2v e» w » (A.31) 

7T 



we rewrite the solution ( |5.1|) in the form 



.(3 cot f + cot v 

z = i—v + 1 in 



7T COt I! — COt V 

cot z f 



2 ^ (C4-C)cott; oo + 2 fC42 



(A.32) 



(C4-C) tan^-fC42 



Putting v = v% in ( |A.32| ), we get 



/3 tanh b^ — i cot i> 

a = 02 + « m 



7r tanh 5 2 + i cot u 



9 1 — TT" COt V 

tanh 2 5 2 - 27r 



1 + tan v 

271" c 



and finally 



cos(a + — #2) = cos2v oo — ^- sin 21?^ ~ cos2t> c 

7T 27T 
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which allows us to evaluate a 

o»2« oc -^ ) S 2 « In f — sin 2^ J . (A.33) 

Appendix B. Direct scaling analysis of the equations on parameters of the large 
N solution 

Six conditions on the length of the cut, a, and the 5 parameters of the map z h- > £(z) 
(which we denote here by a\, . . . ,a n instead of £i, . . . , ( n ) are 



(1) a= dt 

J a 



" t-a 2 

a, \Y(t)\ 



(2) 



£ dt wi\ ( resp - * = C* 



t — a 2 



( 3 ) ^=r dt wm 

Ja 3 \Y(t)\ (B.l) 

(4) a\ + a 3 + a 4 + a 5 = 2a 2 

(5) lua 5 = £^(^-i) (a 5 >0), 

47r • u^/on r 5 ^c(c-a 2 ) 



47r 

(6) — sinh(/3/2) = y rfC 



1^(01 



with Y(t) as in (A3) (where the ai were denoted by Q). While the first three conditions 
are implied by the geometry of the map £, conditions (4) and (5) follow when comparing 
the (from (4.5)) known asymptotics of (, e.g. for z = it + iy (y — > oo; cp. (A4)), to the 
one implied by the integral representation (A2) which says that 



y = inC+ + {^J dt ~ i) ~ lnas ) ~ £~ Q( ai + a 3 + a 4 + a 5 ) - a 2^) - ^ 

n= (B-2) 

with 6 2 = i^( a i + a l + a 4 + a i — n 2 =£i^j=£2 a i a j) already determining d6 p{9) cos(6>). 
With 

, . / Q>41 , V a 54«31 ,/ ^ V a 43 a 51 

aij=a,i — aj, 7 = V a 53 a 4i5 z/ = arcsinw /c = , fe = yl - k z = - 

V a 5i 7 7 

(B.3) 

the conditions (2), (3), (5) and (6) read 





fa 3 i 


a 4 2 






f-° 31 , 


«52 


V a 53 




a32 


\«53 



(B.4) 



7T7 

2a 32 
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(3) -f3 = (^,k')~ ^K(k') 

2 7 ^ \ 7 (B.5) 

2a "n(m,iA + ^ K{k>) 

\ a 53 / 7 



7 

(5) ^±n L fc') = In ( + (B.6) 

7 V a 53 / 7 V 4 / 

(6) fsinh(/y2) = l gW - + "'°^ W . (B.7) 

The scaling limit can be obtained by direct analysis of the equations (|B . 1|) in the limit 
(3 — > 0, a — > ao > 0. Letting 

a 43 = it, a 32 = u, a 53 = w, 

u = , v = , w = , 

^31 &31 a3i 



one has 



(3) .? = 2/ 1 rfs- ? = V + US „ (B.8) 

ys(l — s)(iy — us)(l + us) 



o 



2, /• 1+J. 



tt31 

1 



(1) a = v I ds 1 S (B.9) 

Jo ^s(s + ^)(w + vs)(l-vs) 

1-8 



ds 



wa 3lJo a / s(s+ h )( i + ^)(i__h_ s ) 



7T = — — — I ds- 



V^Ju/w y/s(l-s)( S -u/w)(l + Sw) 

2« r/M r f 1 ds 
~ —?^K{r) + Vto / = 
</o y(l-s)(l + sw) 



where r 2 = (1 + tt) 1 (1 — — ) « 1 — w 



w + l - 



In order to have (3 — > and a finite, we must have u, v — >• 
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■0. (B.ll) 

If one wants to keep, according to ( |B.9| ), a finite in this limit, u/v must go to zero such 
that 

11 

-i\n- = 2a (B.12) 

v 
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finite, i.e. u — > exponentially faster than v/yib (and the ~ sign in (BIO) and thereafter, 
means that such terms are dropped). 

One also finds that — — > (even if w — > 0) as if not, the r.h.s. of ( P.10|) would go to 
zero. So j3 — > 0, a — > ciq > implies 

u,v^0, -->0, — ->0 (B.13) 

together with ([BHD , (|BTT2| ), and 

/3 « eyr. (B.14) 

In order to extract more quantitative information from (|B.10|) consider the equivalent 
condition (|B.4j) , 

(u + v)K(k) = ull ( — k) (B.15) 



where 



1 + u' 



k 2 =^--)(—)< —■ (B.16) 
w) \l + u) 1 + u K ' 



As u — > 0, /c /2 = 1 — k 2 pd — — u — > and we can use some standard expansions for the 
third elliptic integral appearing in ( p.l5| ), e.g. ( 412.01 of [BF]) 

n( 1 A tn' / 1 ; M0 - k) (B.i7) 

1+u 1+-U 



where 



sin6> 



^ w 1+u 



;b.i8) 



and the first terms in the expansion of Heumann's Lambda function Ao(6,k) (904.00 of 
[BF]) are 



A o (0,fc) = - \Ed- -(2K - E)k' 2 ( 6 - sin 6 cos 6) + ... ) « - arcsin J— ^— . (B.19) 

7T \ 4 / 7T Vfi + l 

Inserting ( [H.18| ) into ( |B.15| ) and using 

1 



K(k) « ~2 lnA: + ln4 for &'^0 



one finds: 



1 . f u u \ , \ ir r: L 2 . / to 
In 1 + In4 - vw 1 arcsm 



2 \ W 1 + U J J 2 \ 7T \ 1 + w 

and, using ( P-14j ) 

U , U _2fl(i_2 arcs i n /ZfT) 

I « 16e ^ v * V (B.20) 

to 1 + u y J 
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Apart from converting In it terms into w, 8 dependencies, all other u-dependencies are 
dropped, due to this exponential decay ( [B.21] ). Eq. (B.1.5) can then be stated explicitly 
as an expression for 031 in terms of 8 and w as follows: 

f°° , / 1 1 \ v f°° ds 
In a 5 « / ds I + 



y/(s + a 31 + w)s s + as+wj w J, s ^ (s _ + ^ 

= lim (in (iJ s 2 + sa 5 i + 2s + a 5 i) + In ^ J + — a/u; ( arcsin ( \ — % J + ^ J 
A^oo \ V Jo Ay w \ \l + wj 2 J 

= m 4 + in 1 arcsin H 

a 5 i 2tt V V 1 + */ 2/ 

(B.22) 

Hence 

4 !ri + 2. arcsin (i^5\i 

031 ~ e 4 L w vi+«/J. 115.26) 

1 + w 

The last equation needed to calculate the a{ as functions of 8 — > and A is ( B.7|) , resp. 



4tv sinh(/3/2) , — ; r„/,s 2t><23 rW , s ,^ n .s 

^-L = ^/w(a 31 + u)E(k) + —=^=K{k) (B.24) 

A y/wa 3 i 

as 7 2 = w(a3i + u) and (due to (B.1.4)) 

0405 + aia 3 = —2va 3 . (B.25) 



( B.24 ) can be simplified substantially even without neglecting u-terms, by noting that 
( gig ), (B.17) imply 

' 2r K(k) = VT+ti[l-A (9,k)}. (B.26) 



7t ^^31 

With 



£3 = ^3i (1 _ 2C ^ = A^ (B.27) 

2 4 v h 2tt k ! 



one therefore gets 



2sinh(/3/2) a 3 i , -(2 r- , _ , A 

yiL_L = -£Vl + u -\RbE + (1 - 2v - w - - A ) 

A 4 y 7T / 



a 3 i / 2 r 8 r \ 2 . w 
'—— —Vw + 1 V w — 10 1 arcsm 



4\7T \ 7T J \ n \ w + 1 



(B.28) 



which is an (implicit) equation for w as a function of f3 and A, when inserting ( |B.23| ). For 
tZ; — > it reads 

l-e _/3 8 ,— 

r « 1 - - v 7 ^ - 2u). (B.29) 

A 7T 
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The length of the cut is given by (cp. (B.1.1),(B9)) 



(-/ — (/( i = 7T ( 1 arcsin y -_[_]_ I ' B.30) 



The second of the final scaling eqs. ( |5.6|) follows from (P-28| ) and ( P-30| ) if we neglect 
all terms proportional to (3 or exponentially small terms and use a = Iv^. 

Finally note that the second line of (B.10), via K(r) m | In implies (B.21) 

(shortcutting the argument (B. 15-21)), when using 

i— f 1 ds tv . 1 — w 

yw I — . = arcsin 

Jo - s)(l + sw) 2 1 + w 

and 

11 1 -w 2 . 1 

— I arcsin = — arcsin 



2 TV 1 + W TV sj\ + W 



, 2 . I w 
1 arcsin 



TV V W + 1 



Appendix C. Inverted oscillator: the point (3 = in 

An interesting analytical continuation of our model corresponds to the imaginary 
values of the generator of O e (2) symmetry of the original supersymmetric model. If we 
renormalise e to one it is equivalent to the change /3 — in ( j2.24| ). The corresponding 
saddle point equation reads: 

. f a ,/,/\f n u — u' u — u' + (3 u — v! — (3\ > ~ v 

2A sinu = j-au p[u ) I 2 cot — cot cot 



According to the arguments and results of the paper []14] the inverted twisted matrix 
oscillator describes the compactified c = 1 string, or, in other words, the compactified 
bosonic field coupled to the 2d quantum gravity. So at least the critical regime of c = 1 
string with the typical inverse logarithmic dependence of the physical quantities on the 
cosmo logical coupling should show up at some point. Let us d emo nstrate it in the case 
which we can solve explicitly, namely for (3 = itv. The equation ( |C.1| ) in this case looks as: 



— sin u 



du' . ; v . . (c.2) 



The spectral density is 



2 J_ a sm(u — u') 



p(u) = — V sin 2 a — sin 2 u. 
Itv 



The normalization condition gives 

/ a A f a 2A 

dup(u) = — / duy sin 2 a — sin 2 u = — [E'(sina) — cos 2 a K(sma) 
27T J_ a TV 
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or 

E(k)-k ,2 K(k) = ~, & = sina. (C.3) 
A 

Consider the limit when the eigenvalues occupy almost the whole interval [— 7r, n] allowed 
by the periodicity: a ~ tt, /c 2 ~ 1. In terms of /c' we have the following assymptotics: 

E ~1+ \k' 2 log(4/&/) (C.4) 

K ~ log(4/fc') (C.5) 
By the use of (IClQand (|G3|) we obtain from ( |C^| ): 



7T | log(A - A c ) | 



(C.6) 



for A — > A c = tv. 

For the simplest physical quantity: the derivative of the free energy we obtain: 

f a Xk 2 
F' x = dup(u)cosu= (C.7) 

J -a 2 

from where we obtain the scaling asymptotics typical for the c = 1 noncritical string 



discovered in 23 



4 4 | log( A — A c ) | 

The considered case f3 = in of the c = 1 matrix model corresponds to the Kosterlitz- 
Thouless phase transition point. It would be interesting to study the vicinity of this point 
by generalizing our solution to all imaginary f3. 
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